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Abstract 

We consider the problem of energy-mass transfer from scattering to bound states for a one 
body quantum system subject to the action of a time dependent point interaction in 3-D. Under 
suitable assumptions on the initial state of the particle, we prove a result of local controllability of 
this process. Our proof exploits the finite time asymptotic analysis of fractional integral equations 
and the rank theorem for maps defined on Banach spaces. 

1 Introduction 

A perturbation of the Laplacian supported by a finite set of points in M. d - with d < 3 - defines 

a special case of singular perturbation referred to as point interaction. At a formal level, the associated 
Schrodinger operator can be written as: 

n 

H = -A + Y / ^S{x-y i ) (1) 

i=l 

These operators appeared first in Theoretical Physics during the 30's. They were introduced in 
order to realize a model for the interaction of nucleons at low energy [4]. After, they became a natural 
tool to describe short range forces or "small" obstacles for scattering of waves and particles. 

From the point of view of applications, the main reason of interest of this subject rests upon 
the fact that point interactions often lead to models which are explicitly solvable. It turns out that 
the spectral characteristics (eigenvalues and eigenfunctions) of operators (1), and then all the physical 
relevant quantities related to, can be explicitly computed [5]. This circumstance motivates an increasing 
attention on the application of point interaction models in various sciences, e.g. in physics, chemistry, 
biology, and in technology. 

In this work we build up a point interaction model of a time dependent Schrodinger operator; this 
interaction will be used as a control for the energy transfer between continuous and discrete spectrum 
of a one body quantum system. We will investigate the possibility of finding a time dependence profile 
such that a prescribed part of the energy of a particle, initially placed in a scattering state, moves on 
a bound state in finite time. 

The subject we treat has its natural collocation in the framework of quantum systems control 
theory. Our analysis may find applications in those contexts where short range quantum potentials 
can be used as control tools. 
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Monte S. Angelo 1-80126 Napoli, Italy. Supported through a European Community Marie-Curie Fellowship and in the 
framework of CTS. E-mail: andrea.mantile@dma.unina.it 
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2 Point Interaction Model for Quantum Control Potentials 
and the Main Result 



The rigorous definition of point interactions in Quantum Mechanics - due to F.A. Berezin and L.D. 
Fadeev [3] - rests upon the theory of selfadjoint extensions of symmetric operators. The Hamiltoni- 
ans describing point interactions in the origin of M 3 are defined by the selfadjoint extensions of the 
symmetric operator 

H = -A 

D(H) = C°°(R 3 \{0}) 

In [5] it has been shown that these extensions, denoted in the following with H a , are parametrized by 
a real a 1 ; for any fixed A e M+, the operator H a can be represented as 2 



(2) 



f D{H a ) = [i> G i 2 (M 3 ) | V = 4> X + qG x ; 4> x G ff 2 (M 3 ); A (O) = q(a + } 
\ (H a + X)i/j = (—A + A) cf) x 

where G x is the Green function of (— A + A) 

p -VX\x\ 

G x (x) = -r^r (3) 

47T |^| 

As it follows from (2), every function in D (H a ) is composed of a regular part 4> x G H 2 plus a singular 
term G x . It is worthwhile to notice that this representation is not unique, but depends on the choice 
of A. In particular, for a < 0, we can fix: A = (47r|o;|) 2 , obtaining a null boundary condition: 
</ )A (0) = q( a + 4-7) = 0. In this case the operator domain writes as 



D(H a ) = ^EL 2 (R 3 ) 
where ip a is the normalized Green function 



^^+-^^,^ff 2 (K 3 U(0) = 0,geC (4) 

471V 2 H ) 



72~He- 4 "l"l N 
= ij ( 5 ) 

This representation, holding for negative values of the parameter a, will be extensively used in this 
work. 

We recall the spectral properties of H a : the absolute continuous spectrum coincides, for any a£l, 
with the set (0, +00); the point spectrum is empty for a > 0, or it contains a single eigenvalue: 
X a = —16ir 2 a 2 , for a < 0. The related eigenstate is ip a , defined in (5). 

A time dependent point interaction Hamiltonian H a ^ is defined by fixing a real valued function, 
a(t), which characterizes the time dependence profile. The Schrodinger equation 

ift^ = H a( t)ip 

i){x,$) = Mx) (6) 
1>(x,t) G D(H a(t) ) 



^^From the physical point of view, the parameter a is linked to the inverse scattering length of H a , [5]. 
2 For the definitions and the properties concerning point interactions operators, we refer to the book [5] 
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describes the quantum dynamics generated by H a ^y In [6], D.R. Yafaev has provided an explicit 
expression for the time propagator associated to H a ^ (equation (2.6) in the cited paper) 

i>(t,x) = U t ipo(x) + i Slu(t- s,x)q(s)ds 

(7) 

^ q(t) + AV^Jo ^0fds = iVrif* ^^ds 
with 

UMx)= [ U(t,y-x)^(y)dy; U {t,x) = —^e i ^ (8) 
•Is? (Airit) 2 

Given a € Lf^ c (0, oo), the map (7), acting on the initial state tpo, is unitary in L 2 (M 3 ) and defines, in 
the weak sense, the solution of (6) at time t (sec [6] and references therein). The auxiliary variable q(t), 
solution of the second equation in (7), is usually referred to as the charge associated to the particle. 
It expresses the coefficient of the singular part of the state tl>(x,t). We refer to the Appendix A for a 
study of the charge equation. 

Let us fix, from now on, a < 0. We assume, at t = 0, the particle to be placed in a scattering state 
of the Hamiltonian H & 

^ = ^i}(ff fi ):(^i a ) i2 =0 (9) 
Our purpose is to find a suitable a(t) in the space of control functions 3 

o6flo(0,T) (10) 

such that the control interaction, H a ^ +a , is able to steer the system (7) from ip to a state, i(>(T), 
whose projection (?/>(T), tpa) has a previously fixed value in a neighborhood of the origin of C. 
Replacing a with a + a in (7) and projecting along ip a , we get a complex valued map 



F(a) := (u T <P + l £u(T-s,-)V(a)( S )d S , V>a) i2(R3) 
V(a) = q(t) : q(t) + iVri J* fa(8 ^gL g(a) ds = A^l J* ^§ds 
For a = 0, the system evolves under the action of H a ; in this case we have 

F(0) - {e~ iTH ^, Va) = e lTAa (0, ^) = 
The condition for the solvability of our problem results to be the local surjectivity of the map 

f F{a) = zeC 



(11) 



(12) 



around the point a — 0; i.e. we aim to prove that there exists a neighborhood of the origin in C, Jo, 
and a C 1 -class map g : Io — > Hq(0, T) such that 

F(g(z))=z Vze/o (13) 

The main goal of this work is to demonstrate the following result 



3 The boundary conditions on a being chosen in order to guarantee that at the inital and final times the Hamiltonian 
is H a . 
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Theorem 1 Let <f> be a scattering state of the Hamiltonian H & fulfilling the condition: 

e {7 G S(R 3 ) I ( 7 , ^ a ) L2 = 0; 7 = 7(|£|); 7(0) = } (14) 

being »S(R 3 ) the space of functions of rapid decrease (see e.g. [7]). Then the functional F : Hq(0,T) — > 
C, defined by (11), (10), is a locally surjective map around the point a = 0. 

Remark 2 The inverse problem related to (11), i.e. find a control a such that F{a) has a fixed value 
in C, is connected with the investigation of the energy exchanges in non autonomous quantum systems. 
Although physical intuition suggests that the time dependence of the Hamiltonian can force energy ex- 
changes between the point and the continuous spectrum, the dynamics describing these energy transfers 
is, in general, rather complex (see, for instance, the study of the ionization problem driven by time 
periodic delta interactions in [l]-[2]); in particular, we would like to stress that the local controllability 
of these phenomena is not a simple consequence of the time dependence of the Hamiltonian, but it 
depends from the specific structure of the interaction considered. 

We will adopt here a standard procedure in the analysis of nonlinear systems. First we study the 
regularity properties of the map F : Hq(0,T) — > C and prove that F <E C 1 (Lemma 5). Then we 
investigate the surjectivity of the linearized map d F. To this aim we will study a non controllability 
condition for the linearized system (Section 4); it will be shown that, under the hypothesis (14) on the 
initial state, this condition is never satisfied, obtaining, in this way, a controllability result (sections 
5-7). Then we conclude using a Rank Theorem for functionals defined on Banach spaces (Section 8). 

3 The Nonlinear System 

In this section we investigate the differentiability properties of the functional F(a). We shall denote 
with X and Y two Banach spaces, U an open subset of X, d a F and d^F respectively the Frechet and 
the Gateaux derivatives of the map F : U —* Y evaluated in the point a. A differentiable functional 
F : U — > Y is said to be of class C 1 if the map: 

F' :U -» L(X, Y), F'{a) = d a F (15) 

is continuous. Next we recall a standard result in the theory of differential calculus in Banach spaces 
(Theorem 1.9 in [11]): 

Theorem 3 Suppose F : U Y is Gateaux- differentiable in U. If the map: 

F' G :U ^L{X,Y), F' G {a)^d^F 

is continuous at a* S U, then F is Frechet- differentiable at a* and its Frechet derivative evaluated in 
a* results: 

d a ,F = d%F 

Consider the map V : H^(0, T) -» C(0, T) defined as follows 

[V(a)} (t) - q(t) : q(t) + A^JVi f Ks) +^L^ ds = f(t) (16) 

Jo \/t-S 
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with / e C(0, T). From the estimate (62), in Lemma 10, we have the following bound for the C(0, T)- 
norm of [V(a)\ {t) 

\\V(a)\\ 

C(0,T) < 11/11 C(0,T) 

T(a,T) (17) 
where r(a, T) is a positive finite constant depending on a and T. 

Lemma 4 Let f e C(0,T). TTie map V : i^O,! 1 ) -> C(0,T) de/med &y (16) is of class C 1 . 

Proof. First we prove that V is continuous; let a, (i £ Hq(0, T) and consider the difference V(a)(t) — 
V((3){t): it solves the equation 



q(t) + Wiri 



y(s) + a] q(s] 



which is of type (16); from the estimate (17) we obtain 
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Jo 



\\V(a)-V(P)\\c { o,T) 
moreover, the Sobolev inequality 



< 



ITO)llc(o,T)ll/3-a|lc ( o,T) ^T(a,T) 



\\9\\c(o,t) - C \\9\\m(o,T) 



implies 



\\V(a)-V(P)\\ c{0tT) <C AVni \\V(0)\\cq,T)\\0 ~ <*\\m(o,T)W*r(a,T) 



(18) 

(19) 
(20) 



From (65), it follows that T(a,T) is a continuous function: C(0,T) -> R+. Therefore T(a,T) is 
uniformly bounded for a close to (3 in f/^-norm. Previous remark and the estimate (20) imply the 
continuity of the functional V from H%(0, T) to C(0, T). 

Next we introduce the Gateaux derivative of V in the point a 



d G a V(u) = lim 



V(a + eu) - V (a) 



this limit being in the sense of C(0, T)-topology. The explicit expression of d^V{u) is 4 



d%V{u) = q : q{t) + l\firi 



f 

Jo 



[a(s) + a] q(s) 



ds + 4v7ri 



I 



1 [V{a)\ (s)u(s) 



ds = 



\Jt - s J \/t — s 

This is a linear map : (0, T) — > C(0, T). Making use once more of the estimate (62), we have 



(21) 



\d G a V{u)\ 



C(0,T) 



< 



lin«)llc(o,r)ll«llc(o,T)r(a,T)2T' 



In order to prove the statement of the Lemma we need to show that V : Hq (0, T) — > L(Hy (0, T), C(0, T)), 
defined in (15), is continuous. From Theorem (3), it is sufficient to prove that the map 

V' G : T) - L(Ht(0,T),C(p,T)), t£(a) = d G a V 



'This formula can be obtained by exploiting the result of Lemma 10. 
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dgV-ttjfV 



is continuous, i.e.: linic^ 

Set d^V(u) — d^V(u) = q; from the explicit formula (21), it follows 

, s ./—■[* [a(s)+a] q(s) , i 7 I' 
Jo y/t-s Jo 



0, where ||-|| is the operator norm in L{Hq (0, T), C(0, T)). 

f 

* {[V((3)\ (s) - [V(a)\ (a)} u(s) [a(s) /?(*)] [V{(3)\ (a) 



This equation is still of the type (16); from the estimates (17) and (19) we get 



-ds 



\l\\c(0,T) 



< 



4V7ri 



then the operator norm 



"llffW) WW) - V{a)\\ c{0 T) + \\a- P\\ m(0 , T) ||/|| c( o,d r(/J,T) T(a,T) 
d^V — dpV can be bounded as follows 



\d°V-dfjjV\ 



sup \\d^V(u)-d^V(u)\\ c{QT) < 



llu|| = l 



< 



4V7ri 



\\V{(3) - V{a)\\ c(0tT) + \\a- /3|| ffl(0 , T) ||/|| C(0 , T) T(/3,T) T(a, T) 



As already noticed, F(a,T) is uniformly bounded when a is close to j3 in iJ 1 (0, T)-norm. Then, the 
continuity of the map V allows us to conclude that 

lim |||d?V-dfV||| =0 

■ 

Our next task is the study of the regularity properties of the map defined in (11) whose explicit 
expression is recalled here: 



F(a) := (u T <t>{x)+i U(T - s,x)V(a)(s)ds, ifefe)) 
V{a) = q(t) : q(t) + ^ J* MAjgjlsl ds = 4^/* ^ 
Lemma 5 The map F : H%(0, T) -> C defined by is of c/ass C 1 . 



(22) 



Proof. First we show that F is a continuous map. By inverting space and time integrals in (22), the 
scalar product (j, U(T — s,x)V(a)(s)ds, ips(x)j ^ ^ is given by 



if U(T-s,x)V(a)(s)ds,if> a (x)) =if f U(T - s,x)iP & {x)d 3 xV(a){s)ds = 

Jo I ,,,„,> Jo ii 3 



L 2 (R 3 ) 



= i / L/- T -Wa(Q)V(a)(s)ds 
Jo 

Let a,(3 € -ffo(0, T); making use of the previous expression, the difference F(a) — F{(3) can be written 
as 

F(a) - F(/3) = i / U T _MQ) (V(a)(s) - V(J3)(s)) ds 
Jo 
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Therefore, it satisfies the estimate 



\F(a) F(P)\ < \\V(a)(s) - V((3)(s)\\ c ^ T) £ Pt-.^M *s (23) 

The continuity of F then follows directly from Lemma 4. Moreover, from Theorem 3 and definition 
(15) it follows that F is of C 1 class if the map 



is continuous, i.e.: lim^-^ 
a and acting on u is 



F' G : Hl(0,T) - L(^(0,T),C), F' G {a) = d G F 

= 0. The Gateaux derivative of F evaluated in the point 



d G F-d G F 



d G F{u) = i (j U(T - s,x) d a V(u)(s)ds, Vafe) 



(24) 



By inverting space and time integrals, we obtain 



d G F{u) = i / Ut-sMQ.) d G V(u)(s)ds 
Jo 



(25) 



Consider the difference d G F(u) — d G F{u): from (25) it results 



\d G F(u)-d G F(u)\ = 



f Ut-sMQ.) (d G V(u)(s) - d G V(u)(s)) ds 
Jo 

< \\d G v(u) d G v(u)\\ c{0T) £ \U T -sMQ)\ d S 



and 



\\\d G F-d G F\\\= sup \d G F(u)-d G F(u)\<\\\d G V-d G V\\\ [ \U T -sMQ)\ds 

ueH*(0,T) Jo 

IMI=i 

The continuity of the map F' G easily follows from the continuity of the map V' G . ■ 



4 The Linearized System 

We consider the map: 



doF(u) = z 

ueff 1 (o,T) 



(26) 



where d^F is the Frechet derivative of functional F evaluated in a = 0. Our aim is to prove the 
following result: 

Theorem 6 Under the assumptions of Theorem 1, the map defined by (26) is surjective. 
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The proof of Theorem 6 will be given in Sections 4-7 following an ad absurdum procedure. 
First we write the functional (26) into an explicit form: 



' d Q F(u) = i (/ T U(T - s, •) d V(u)(s)ds, Va 



L 2 (R3) 

(27) 

d V(u) = q{t) : q(t) + AV^aJ* $Lds = -*Vri J* ^Bisl ds 

Here the dependence of doF(u) from V(0)(t) may be emphasized by making use of the relation (77) 
(in Appendix A), plus the Fubini Theorem (in order to exchange time and space integrations in (27)); 
proceeding in this way, we get 



d F{u) = i f Ut-sMQ) d V{u){s)ds = 
Jo 

= Ani'i \ U T - s ipa{0) / G(t- s')V(0)(s')u(s')ds'ds 
Jo Jo 

Applying the Dirichlct formula to the double integral we finally obtain 

duF{u)=Ami V(0)(s')u(s r ) U T - s ^a(0) G(s - s')ds ds' (28) 

JO is' 

Let us suppose the map d F : Hq(0,T) — > C to be non-surjective; then d F(u) should have 
a constant direction in the complex plane for any u e Hq(0,T). We can express this as a non- 
controllability condition 

3 ° e ° : { c ^ v(p)W XJ UtsMQ) g(s sides' = o yue H ^ T) 

where '•' indicates the scalar product in C. In particular, being this condition true for any real valued 
u e Co°(0,T), the fundamental lemma of calculus of variations implies 

3 ° G C : { C % m f? U T . sM 0) G (s - t)ds = V< ^ [ °' T] (29) 
Making use of (80) - in Appendix A - we get 

5 Small Time Asymptotic for the Charge 

In this Section we consider the asymptotic behavior of the charge V(0) for t — ► 0. We shall denote 
with o(t) and 0(t) two complex valued functions of the real variable t satisfying the conditions 

lim^ = (31) 
t-o t 
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limsup ^ < oo (32) 
t->o t 

According to our hypothesis on the initial state, (14), and to definition (11), the function V(0) is 
the solution of the equation 

q(t) + AV^a f -^ds = 4V™ f ^Mds 

Jo V* ~S Jo y/t - S 

therefore, thanks to (76), it may be represented in the form 

V(0)(t) = 4ttV7 f G(t - s)C/ s7 (0) ds (33) 
Jo 

Its small time behavior is connected to the limiting behavior of G(t) and U t j(0) for t — > 0. In order 
to study this problem we need the following Lemmas 

Lemma 7 Let 7 = 7(|a;|) be a radial function belonging to the space of functions of rapid decrease 
S(M. 3 ) (see e.g. [7]). If we assume that 



* (34) 



f d n 
Do = |n e N : ^7® 

is a non empty set, then the function U t "f(0) admits the power expansion: 

U a (0) = a m t m + 0(t m+1 ); a m ? (35) 
An (-i) m 



/ k 2m+2 Tj(k) dk (36) 
Jo 



( 27r )2 mi 

where T denotes the Fourier transform, m — minDo o,nd 0(t) € C°°[0, +00). 

Proof. The Fourier transform operator, T, is an homeomorphism of the space 5(M 3 ) in itself. It acts 
on Ut7(Q) as follows: 

TU tl {k) = T^k) e~ ikH 
Then, using J 7-1 , we can represent E/t 7 (0) in the form 

a r+co 

U a (0) = r / k 2 T-i{k) e- lkH dk (37) 



(2tt) 2 Jo 

From the regularity assumptions on 7, we have: T^(k) <E S^R 3 ) and C/*7(0) € C°°[0, +00). Setting 
m = minDo, the Taylor's expansion of 11^(0) up to order m in a right neighborhood of the origin is 

U a (0) = a m t m + O(t m+1 ) 

with 0{t) e C°°[0, +00). The coefficient a m is different from zero - due to the hypothesis (34) - and 
explicitly given by ^^U t j(0)\ t=0 ; form (37), this quantity is 

4 (A\ m r+00 
^ ( l) 1 k 2m+2 ^(k)dk 



( 27r )l m! 
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Lemma 8 Let G{t) be given by (75); then it admits the following representation 

G(t) = -L+b + O(t); b = -4naVi (38) 
virt 



where 0{t) e C°°[0,+oo). 

Proof. The proof easily follows from the analytic properties of the erfcit) function ([8], relation 7.1.5, 
pag 297 ) ■ 

We will use these results to get an expansion in power of t? for the charge (33). If we assume 
condition (34) to hold, Lemma 7 may be applied to our case, with the only restriction m ^ due to 
the fact that the boundary condition 7(0) = implies Ut"f(0) — for t = 0. By substitution of (35) 
and (38) into (33) we obtain 

rt s m r-t rt Qf s "i+1\ 

V(0)(t) = AVwiam / ds + iirVi a m b / s m ds + 4Viri V ds+ 

Jo sjt - s Jo Jo Vt - s 

+ 4irVib f 0{s m+1 )ds + 4irVia m f 0(t - s) s m ds + AirVi f 0(t - s) 0(s m+1 )ds 
Jo Jo Jo 

An explicit calculation of the first terms in this expression leads us to the following expansion 

V(0)(t) = A m a m Vit m+ 2 + Bm a m V~ib t m+1 + 4V™ / \ ' ds+ 

Jo Vt - s 

+ 4TrVib [ 0(s m+1 )ds + 4TrVia m [ 0(t - s) s m ds + iirVi [ 0(t - s) 0(s m+1 )ds 
Jo Jo Jo 

where A m and B m are strictly positive real constants. By definition (32), it exists C such that 
1 0(s fc ) I < C\s k \ < C\t k \ for any |s| < \t\; using this estimate, the following relations are easily 
obtained 

/ 0(s m+1 )ds = 0(t m+2 ) (39) 
Jo 



f 0(t-s)s m ds = 0(t m+2 ) 
Jo 



(40) 



0(t - s) 0{s m+1 )ds = o{t m+2 ) (41) 

* 0(s m+1 ) -h 

>ds = 0{t m+ i) (42) 



lo Vt - s 

and the small time asymptotic representation for V(0) can be rewritten as 

V(0)(t) = A m a m Vit m +i + B m a m V~ibot m+1 +0(t m+ l) (43) 
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6 The Non Controllability Condition in the limit t — > 

Here we study condition (30) in a neighborhood [0, 5) of the origin with 6 <T. 
Let us first set the condition (30) in the equivalent form 

3Ke [0,2tt) : arg (V(0)(t)e- l(T - t)A5 ) =K Vie[0,T] (44) 
where argz is defined modulus 2ir. Making use of (33), we see that 
arg(V(0)(t)e- 4(T -* )A *) = 

= arg (^TtViJ^ G(t - s)C/ s7 (0) ds e -<(r-t)A« j 
Then, condition (44) implies 

3 K e [0, 2vr) : arg (^iryfi J G(t - s)U s j(Q) ds e -i(r-t)A a ^ =K Vie [0, T] (45) 
In order to analyze (45), we need the following Lemma: 

Lemma 9 Let a m and bo be defined by (36) and (38) respectively. Under the assumptions of Lemma 

7 and Lemma 8, the function 



arg (Wi G{t - s)U sl {Q)dse-^ T - t ^ 



admits the small time expansion 

rt 



arg^4?iVijf G{t- s)U s -/{0)dse- l{T ^ x ^ = 

= arg (a m Vij + c sin(arg b )t* - (T - i) A a + o(t* ) (46) 

wzi/i c = ^ | & 1 ■ 

Proof. Set z\ — pie'^and z 2 = P2e lip2 ; the first order Taylor expansion of arg(zi + z 2 ) w.r.t. the 
ratio e = ^ about the point e = is 

arg(zi + z 2 ) = <pi + sm(f 2 - <pi) e + o(e) (47) 
Using (33) and (43)we have 

arg^Wijf G(t - s)U s -y(0) ds e^-^*^ = 

= arg (A ro a TO V~it m+ i + B m a m V~ib t m+1 + 0(t m +§)) - (T - t)A a (48) 
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Using (twice) relation (47), the right hand side of (48) can be expanded as 

arg (A m a m V~it m+ i + B m a m V~ib t m+1 + 0(t m+ i)) = 
= arg (^A m a m Vit m+ i + B m a m V~ib t m+1 ^ + O(t) = 

= &vg(a m Vtj + csin(argfeo) t% + o(t%) (49) 

Equation (46) is a straightforward consequence of (48) and (49). ■ 

Lemma 9 leads us to an asymptotic formulation of the non controllability condition for small time. 
From relations (45) and (46), indeed, we have 

K = arg (a m Vtj + -jj= i* - (T - t) X & + o(t' ) (50) 

where the explicit value bo = —4naVi has been taken into account. Recalling that c ^ 0, relation (50) 
is an evident contradiction we obtained supposing the system (26) to be non-surjective. This concludes 
the proof of Theorem 6 for all choices of initial states satisfying condition (34) of Lemma 7. 

In the next section we will study an extension of the proof to those cases in which Lemma 7 does 
not applies. 



7 Finite Time Asymptotic for the Charge and Proof of 
Theorem 6 

If condition (34) does not hold, we may still recover our results by changing the point in which we 
perform the expansions of expressions (33) and (44). 

To this concern, we consider a radial function 7 = 7(|a;|) in the space S^IR 3 ). From (37) in Lemma 
7, it follows that, Ut"f(0) is a C°°-class function represented by 

Ua(Q) = t k 2 Tj(k)e- tk2t dk (51) 

(2tt) 2 Jo 

where T denotes the Fourier transform in i 2 (M 3 ). Let us define / e L 2 (— 00, +00) 

f(y) = { {2i: ^ h V*Fl{v h ) V>° (52) 

Setting k 2 = y in the integral (51), we can express U t j(0) as the Fourier transform of / 

r+00 

Ua(Q) = f(t) = / f(y)e-^dy 
Jo 

Making use of relation (53), it is possible to extend £/*7(0) to the complex plane as follows 



(53) 



r+00 

UMO) - / l{y)e~ tyz dy- z = t + is (54) 
Jo 
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It is well known that, for any / G L 2 (0,+oo), equation (54) defines an holomorphic function in the 
lower complex half plane s < ([9], Section 19.1). In order to study the limit of £47(0) as z approaches 
the real axis, we notice that this function can be expressed as the Fourier transform of the product of 
f(y) times e^ s e L 2 (-oo,+oo) 

tfz7(Q) = (f(v)-e\y\'){t) 

therefore we have 

U t+is l(0) = f * (eH«)(t) = - / U t - tn (0) - ' 1 dt' (55) 

J-oo \t'\ + S A 

with 



l - 2 = (^ s )(t) (56) 



* \t'\ +s 

Taking into account the continuity of L/t 7(0), it can be shown that £47(0) is continuous in the set 
s < 0. 

So far we obtained that Ut+i S "f(Q) is holomorphic in {t € M; s < 0} and continuous in the closure 
{ieK; s < 0}. For 7 ^ 0, this implies that the zeroes of function U t j(Q) have to be isolated points. 
If 7 belongs to the domain (14), the origin of the real axis is a zero of £/*7(0); therefore a time t > 0, 
arbitrarily near to the origin, exists such that: £^ 7(0) 7^ 0. Proceeding as in Lemma 7, it is possible 
to obtain for the function U t j(0) the following power expansion around t 

£47(0) = a + 0((t - t )); a ^ (57) 

with a = £4,7(0) and 0(t) e C°°[0,+oo). 

Next, we observe that a simple change of variable 



r = t — to 

and the use of (76), provide us with an equation for the charge when the initial time t = to is assigned 

V(0)(t) = AirVi f G(t - s) U to+sl {0)ds (58) 
Jo 

Using (38) and (57), we get the power expansion of the of the non controllability condition in a right 
neighbourhood of the point t = to 

K = arg (a Vtj + -jj= (t - t ) h - (T - t)X & + o((t - to)*) (59) 

with c ^ 0. As in the previous case, this relation constitutes a contradiction obtained supposing the 
system (26) to be non-surjective. 

This concludes the proof of Theorem 6. 



8 Proof of the Main Result and Final Remarks 

So far, we succeeded in proving that the functional F(a), defined by (11)-(14), belongs to C^if^O, T), C) 
and its derivative, evaluated in a = 0, is surjective. The Rank Theorem (see e.g. in [10], page 336 
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Theorem 34), then, applied to our case, implies the existence of a neighborhood of the origin in C, Jo, 
and a C 1 -class map g : I — > Hq(0, T) such that: 



F(g(z))=z Vze/o 



(60) 



This concludes the proof of Theorem 1 . 

Our main remark is about the assumptions (14) on the initial state. By considering functions with 
radial symmetry, we are taking into account only those scattering states which have a null projection 
along all spherical harmonics excepting the first one. In this choice there is no loss of generality. 
Indeed, those scattering functions whose expansion in spherical harmonics is: 



Then it follows from definition (7), and the uniqueness of solution of the charge equation, that a 
particle, initially placed in such a scattering state, <j>, and subject to the action of any Hamiltonian 
of type H a i t \ , results to have a null charge and to evolve under the action of the free propagator; in 
other words, starting from this initial condition, the particle doesn't feel the interaction at all. In this 
case any transfer of energy is physically impossible. On the other hand our model can be applied to a 
realistic situation in which an incoming particle - described by a wave function of type: ip — 4>(r)f($, ip) 
with /(■$, ip) null outside a cone - is partially trapped into the attractive potential described by H a . 

In conclusion, we have proved the local controllability of a process of energy-mass transfer, from 
scattering to bound states, for a one body quantum system under the action of a time dependent point 
interaction. 

Further development of this studies may concern the global controllability of the same process, as 
well as the inverse problem of finite time ionization. 
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exhibit the following characterization: 



U t (j)(0) =0Vt 
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Appendix A: The Integral Equation 



We recall some basic properties of the integral equations we deal within this work; some of the relations 
already used in the previous sections will be here obtained. A detailed analysis of fractional integral 
equations, which arise in the framework of time dependent point interactions in Quantum Mechanics, 
is given in [1]. 



Lemma 10 Let a, f e C(0, T); the equation 

rt [a{s) + a] q(s) 



q(t) + 4 



ds = f(t); a € 



10 y/t-S 
has an unique continuous solution such that: 

IMIc(0,T) < ll/llc(0,T) r ( a > T ) 

where T(a, T) is a finite positive constant depending on a and T. 

Proof. The solution q(t) may be formally expressed using the Picard series: 

9(*) = E^9n(t) 

qo(t)=f(t) 

qn (t)=4V™Jo la{s)+ ^ a - l{s) ds 
which admit the following estimate: 

+oo 

A ' < " f\\C(0,T) 



(61) 



(62) 



(63) 



+oo 



i+Eh 



n=l 



l a + a llc(0,T) 



(64) 



with: 



n odd 



r(i) nl! 

being nil = 1 • 3 • 5 • • • n for n odd. The strong infinitesimal character of the sequence A n allow this 
sum to converge uniformly for any a £ C(0, T) and for any finite time interval [0, T\. The existence of 
a unique solution of (61) satisfying the estimate (62) directly follows from (64), with 



T(a,T) = 



i + E* 



HI 



n=l 



l« + «llc(0,T) A n ^T^ 



(65) 



In connection with the charge equation in (7), we claim the following result: 

Lemma 11 Let ip e D(H a ) where H a is the Hamiltonian associated to a 3-D point interaction placed 
in the origin and a G M; then, the function: 



o \Jt - s 



ds 



(66) 



is continuous. 
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Proof. First we consider the case a < 0. 

From definition (4), any function ip € D(H a ) is the sum of a regular part plus a bound state term: 



ipo = <P + 



471-^2]^ 

Then the function (66) disparts in two contributions: 



i> a , ip e H 2 (W), ip(0) = 0, q e C 



f u -4M d8 =f 

Jo V* - s Jq 



U s p(0) 



ds 



Vt-s ^s/2\ 



= / 

2 la 7o 



y/t^S~ 



ds 



(67) 



(68) 



We want to prove that (68) defines a continuous function on finite time intervals. To this aim we 
consider the two contributions of (68) separately. 

The first term of the second member is the one half integral of Utp(0). Using for the state <p a 
representation in terms of spherical harmonics, holding for all L 2 (R 3 ) functions, we have: 

+oo / 
1=1 rn— — l 

From the orthogonality relations for YJ m , we have 

umo) = u tX (o) 



X € H 2 (R 3 ) denoting the radial part of p. Thus, the function U t ip evaluated in x = may be expressed 
by the following Fourier integral 



UMQ) 



1 



(2tt) 2 JR 3 



Ik- 1 t 



dk 



(69) 



The Fourier transform of the radial function x € H 2 (R 3 ) has the following characterization: k 2 Tx{k) € 
L 2 (K 3 ). Using this relation and the Schwartz inequality, we get an estimate for the Z^-norm of Tx 



ll^xlli= / 

Jm 3 



\f X (k)\(l + k 2 ) 
1 + k 2 



d 3 k< + k 2 ) T x (k)\ 



1 



1 + k 2 



from which it follows that Tx € L 1 (R 3 ). This result guarantees that the function (69), as well as the 
first source term f Q * ds, are continuous. 

The second source term of (68) may be evaluated explicitly by using the Laplace transform operator 
£; from definition (5), a direct calculation shows that 



Lo 



V~t~s~ 



ds 



(P) 



2tt lal 



P 



Here we recall that 



i pi — Airay/i 



(70) 



(71) 
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Then the second contribution of (68) is 



ds 



47r^2 \a\ Jo Vt- s A\fni 



=e lW7T2a2t erfc(4ir\a\Vii) 



(72) 



which is a continuous function bounded in M + . 

The same result holds in the case a > 0, where the only difference consists in the fact that the 
operator domain does not include any bound state of H a . ■ 

Relation (62) allows us to obtain an estimate for the solution of charge equation in (11) 



l9llc(o,T) 



< 



4V7ri 



o 



ds 



T(a,T) 



(73) 



C(0,T) 



where the boundedness of second member is assured by Lemma 11. 
Solving the charge equation 

In the particular case a = 0, the charge V(0)(t) satisfies a fractional integral equation: 



q(t) + A\firia 

Jo 



g(«) 



ds = 4v iri 



f 

Jo 



ds 



(74) 



whose solution may be explicitly expressed as a functional of the source term; the Laplace transform 
of (74), indeed, gives 



*> I 1 + — J =^F C {Utm> lv) * *> = ^^Twi 



=£(W(fi))(f>) 



and taking into account the relation: 

1 1 



(4nayfi)e i(4 ™ )2t erfc{4iTaVii) = G(t) 



an explicit expression for the charge is obtained: 

V(0)(t) = AirVi f G(t - s) U s <j>(0) ds 
Jo 

Due to its similar structure, an analogous expression holds for equation (27): 

d V(u){t) = AirVi f G{t - s) u(s) V(0)(s) ds 
Jo 



(75) 



(76) 



(77) 



When (j> coincide with the bound state ip a of the point interaction Hamiltonian H a , the charge 
equation (74) is explicitly solvable. Indeed, the time evolution of a quantum particle starting in the 
state ipa and moving under the action of H a , is 
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Comparing this expression with (4), we get the explicit form of the charge 



l(t) = 47T y/2\E\ , 



-itX 5 



Relations (76) and (78) imply 



Vi [ G(t- s) U s ip a (0) ds = y/2\E\ 
Jo 



-it\ B 



(78) 



(79) 



Replacing t with T — t, the previous formula can be written as 

V [ G{T-t-s) U s ip a (0)ds = V^e^ (T -* )Aa 
Jo 

Moreover, setting s' = T — s in the integral at the l.h.s. and dividing by V, we get 

J G(s' -t)U T - s >M0)ds' -- 



2\a\ 



-i(T-t)\ e 



(80) 



This relation has been used in order to obtain condition (30). 

In order to justify the use of Laplace transform in deriving relations (76) and (77), our next task 
is to prove the following 

Lemma 12 Let q{t) be the solution of the charge equation: 



q(t) + iVri / 1 — - J — ds = aVki / 
Jo Vt - s Jo 



o Vt- s 



ds 



with a € L°°(K) ; ipo € D(H a ^). Then, the Laplace transform Lq(p) exists and is analytic at least ir. 
the open half plane ofC defined by the condition: 



peC:Rcp> 167r 2 ||a||^ 
Proof. Consider the function: q'(t) — e~' pt q{t) with p € C and Rep > 0; it satisfies the equation: 



q{t) + A\hH 



L 



* a(s) e -P(*- s ) q(s) , . ,—. f* e~ p[l ~ " : 
Vt^s 



-ds = 4Vri 



1 

Jo 



f(s) ds 



Vt~s~ 

where the function / in the nonhomogeneous term is given by: 

f(t)=e-P t U t M0) 

Using relations (67) and (69) of the previous proof, we may write (83) in the following form: 

1 



f(t) = e-? t 



T**p r (k) e 



ik 2 t 



dk 



(27r) 2 -'k 3 47r^/2 \a 

with !Fip r (k) G L 1 (R 3 ); again an explicit calculation shows that: 



:%(0) 



(81) 

.: 

(82) 
(83) 

(84) 



%(o) - c- 1 



2\a\ 



i (pi + 47r \a\ Vi 



2\a\ 



)= - (4tt \a\ V)e t{A ™ )H erfc^n \a\ Vtt) 
nt 
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from which we deduce that f(t) <G i 1 (0, +00). 
Next we apply the Young's inequality: 

H/*0lli<ll/lliNli 

for convolutions of i 1 (0, +00) functions, to the equation (82), and obtain the following estimate: 



q'W, 1-40F|M| 



+ OO 



-pt 



Vt 



dt < I 



r+00 


e -pt 








dt 


Jo 


Vt 



which provides an effective bound for the norm Wq'^ if the coefficient ^1 — 4y / 7r WaW^ J 
is positive. Recalling that, for Rep > 0, holds the equality: 



+00 



Vt 



dt) 



r+00 


e -pt 








dt= y 


Jo 


Vt 



Rep 



we get the condition: 

1 — 4-y/TT |M 



r+00 


e -pt 








dt^j 


Jo 


V 



dt) > => 1 > AVt 



Rep > Wit \\a\\ 



(85) 



Following the same line, it's easy to prove that, if the condition (85) holds, the partial derivatives of 
the function q' w.r.t the real and the imaginary part of p - both given by: 



te- pt q(t) 



are bounded by integrable functions of t: 

/ \te-"*q{t)\dt< / te~^ ^ +s ) l q{t) 
Jo Jo 



dt < co 



Then e~ pt q(t) is C 1 integrable w.r.t. t G [0, +00) for any p in the domain (81); moreover, in the same 
hypothesis, its partial derivatives w.r.t. p are bounded by measurable functions of t. This allows us to 
conclude that the Laplace integral: 



Cq( P ) 



q(t)e- pt dt 



defines a C 1 class function for p in the domain (81). 
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